Accepted for Publication in AJ 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



HI SELECTED GALAXIES IN THE SLOAN DIGITAL SKY SURVEY II: THE COLORS OF GAS-RICH 

GALAXIES 

Andrew A. West 1 ' 2,3 , Diego A. Garcia- Appadoo 4 ' 5 , Julianne J.Dalcanton 3 , Mike J. Disney 4 , Constance M. 

Rockosi 6 , Zeljko Ivezic 3 

Accepted for Publication in AJ 



o 
o 

(N 
-t— > 
O 

O 

(N 

o 

u 

43 
6 

CO 



> 
d 

ON 
O 



ABSTRACT 

We utilize color information for an Hl-selected sample of 195 galaxies to explore the star formation 
histories and physical conditions that produce the observed colors. We show that the HI selection 
creates a significant offset towards bluer colors that can be explained by enhanced recent bursts of star 
formation. There is also no obvious color bimodality, because the HI selection restricts the sample 
to bluer, actively star forming systems, diminishing the importance of the red sequence. Rising star 
formation rates are still required to explain the colors of galaxies bluer than g — r < 0.3. We also 
demonstrate that the colors of the bluest galaxies in our sample are dominated by emission lines 
and that stellar population synthesis models alone (without emission lines) arc not adequate for 
reproducing many of the galaxy colors. These emission lines produce large changes in the r — i colors 
but leave the g — r color largely unchanged. In addition, we find an increase in the dispersion of galaxy 
colors at low masses that may be the result of a change in the star formation process in low-mass 
galaxies. 

galaxies: photometry — galaxies: fundamental parameters 
- surveys — radio lines: galaxies 



Subject headings: galaxies: evolution 
— galaxies: general 



1. INTRODUCTION 

The star formation history, metallicity, and current 
star formation rate all contribute to the observed colors 
of a galaxy. Remarkably, these factors work in concert 
to yield a fairly well defined locus in galaxy color-color 
space (Strateva et al. 2001). Deviations from this locus 
as well as the morphology of the locus itself, can lead to 
significant insight into the underlying processes taking 
place within galaxies. Many previous studies have uti- 
lized the broad-band colors of galaxies to investigate the 
star formation histories and metallicities of galaxies (e.g. 
Tinsley 1972; Searle et al. 1973; Tinsley & Gunn 1976; 
Balcells & Peletier 1994; Roberts & Haynes 1994; de Jong 
1996; Bell & de Jong 2000; Galaz et al. 2002; Gavazzi et 
al. 2002; Bell et al. 2003; MacArthur et al. 2004; Zackris- 
son, Bergvall & Ostlin 2005; Driver et al. 2006; Skibba 
et al. 2008). 

The advent of large surveys such as the Sloan Digi- 
tal Sky Survey (SDSS; York et al. 2000) and the Two 
Micron All Sky Survey (2MASS; Skrutskie et al. 2006) 
has created a wealth of uniform data and the statis- 
tical foothold to investigate the bimodality of galaxies 
(Baldry et al. 2004; Kauffmann et al. 2003, 2004), lumi- 
nosity function (Blanton et al. 2001; Ball et al. 2006), 
and the average properties of nearby galaxies (Blanton 
et al. 2003a, 2003b, 2005; Geha et al. 2006; Mailer et 
al. 2008; Skibba et al. 2008). While the large optical 
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and infrared surveys have made large contributions to 
our understanding of galaxy evolution, they only trace 
the stellar component of galaxies and do not trace other 
baryonic material such as cold gas. 

Galaxies in the local universe span a range of star 
formation histories - from blue, gas-rich, low-surface- 
brightness (LSB) galaxies that are slowly turning their 
gas into stars with low star forming efficiencies, to red, 
gas-poor galaxies, that have formed the bulk of their 
stars in the past. Stars dominate the visible light output 
of most galaxies, and thus galaxies detected by tradi- 
tional optical or infrared imaging have well developed 
stellar populations. In contrast, the natural way to iden- 
tify gas-rich, less evolved galaxies is by their 21 cm HI 
radio emission. Aside from its importance for global star 
formation, a sample of galaxies with both gaseous and 
stellar information allows for a more complete census of 
the local baryons (a constraint vital to the calibration 
of n-body simulations; Governato et al. 2007; Brooks et 
al. 2009). 

Previous studies have combined large HI surveys with 
optical and infrared samples, namely the Arecibo Duel 
Beam and Slice Surveys with the Two Micron All 
Sky Survey (2MASS; Jarrett et al. 2000; Rosenberg et 
al. 2005) and the merging of HIPASS with SuperCOS- 
MOS (Hambly et al. 2001; Hambly, Irwin & MacGilli- 
vary 2001b; Hambly et al. 2001c; Doyle et al. 2005). 
Rosenberg et al. (2005) were able to probe the baryonic 
content of a large sample of galaxies, but were limited 
by the shallow depth of 2MASS, which does not have 
data for many of the LSB galaxies in the sample. The 
HIPASS/SuperCOSMOS sample of Doyle et al. (2005) 
contains optical data for more than 3600 HI selected 
galaxies but also suffers from the shallow depth of the 
SuperCOSMOS optical data. 

Recent studies have combined the Parkes HI Equa- 
torial Survey (ES) with the SDSS (Disney et al. 2008; 
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Garcia- Appadoo et al. 2009; West et al. 2009; hereafter 
W09). The deep SDSS optical data provides informa- 
tion about the stellar content for all ES HI sources where 
the two surveys overlap. In addition, the uniform, accu- 
rate, and well-calibrated photometry of the Hl-selected 
ES/SDSS sample allows for a more detailed exploration 
of the factors affecting galaxy colors, particularly for 
galaxies with large reservoirs of gas. The optical col- 
ors are reasonably sensitive to age, although IR colors 
are needed to constrain metallicity (Bell et al. 2003). 
There is thus an unavoidable degeneracy between age 
and metallicity when using only the optical colors avail- 
able with SDSS (Bell & de Jong 2001; Bell et al. 2003). 
Our sample does not have a complete set of near IR coun- 
terpart data and some of our results will reflect this lim- 
itation. 

In this paper, we briefly describe the ES/SDSS sample 
in §2 and examine the colors of Hl-selected galaxies by 
comparing them to stellar population synthesis models 
(§3.1) and the colors of optically selected galaxies (§3.2). 
We also investigate how line emission affects the broad- 
band colors of gas-rich galaxies (§3.3). We demonstrate 
an increased dispersion in the colors of galaxies at low 
masses, and investigate its possible origins (§3.4). We 
summarize and discuss our results in §4. 

2. DATA 

The HI data for our sample come from the Parkes 
Equatorial Survey (ES; Garcia- Appadoo et al. 2009), a 
blind HI survey of the southern sky that covers a velocity 
range from -1280 to 12700 km/s with an RMS noise of 13 
mjy, using the multibeam receiver on the 64m radio tele- 
scope in Parkes, Australia. The ES, which is described 
in detail in an accompanying paper (Garcia- Appadoo et 
al. 2009), circles the celestial equator between -6 < S < 
+10 and contains over 1000 sources in 5738 square de- 
grees. The raw data forms part of the HI Parkes All Sky 
Survey (HIPASS; Barnes et al. 2001; Meyer et al. 2004; 
Zwaan et al. 2004; Wong et al. 2006), carried out with the 
same instrument over the entire sky between -90 < S < 
+25. However, the ES fields were searched much earlier 
(Garcia- Appadoo et al. 2009) in readiness for compari- 
son with SDSS data. While the search techniques were 
much the same as those of the HIPASS team and rely 
heavily on their procedures, the source lists are not iden- 
tical. For example, the completeness limit of the ES list 
is 30% fainter than the HIPASS limit. This difference 
is mainly due to our ability to follow-up and confirm a 
higher proportion of the fainter sources, a process that 
would be impractical with the larger survey. The velocity 
resolution of the ES HI spectra is 18.0 km s _1 and the 3tr 
HI mass limit of the survey is 10 6 x Dm pc Mq, assuming 
a 200 km s _1 HI galaxy profile. For detailed descrip- 
tions of the data acquisition, calibration and reliability 
see Garcia- Appadoo et al. (2009) and the HIPASS analy- 
sis contained in Barnes et al. (2001), Meyer et al. (2004), 
and Zwaan et al. (2004) . 

The optical data for this study come from the Sloan 
Digital Sky Survey (SDSS; York et al. 2000; Gunn et 
al. 1998; Fukugita et al. 1996; Hogg et al. 2001; Smith et 
al. 2002; Stoughton et al. 2002; Pier et al. 2003; Ivezic 
et al. 2004; Gunn et all. 2006) Data Release 2 (DR2; 
Abazajian et al. 2004) sky area. The DR2 area is 3324 
deg 2 , about half of which overlaps with the equatorial ES 



region discussed above. Because the SDSS photometric 
software (PHOTO; Lupton et al. 2002) was not optimized 
for angularly large galaxies, all of the photometry pre- 
sented for the ES/SDSS survey has been re-processed us- 
ing the techniques described in W09. The W09 study de- 
scribes the optical sample selection and matching to ES 
and presents the updated SDSS photometry for 195 Hl- 
selected galaxies using both a sky subtraction procedure 
that is optimized for angularly large galaxies and a cor- 
rection for the over-deblending ( "shredding" ) of nearby 
galaxies by the SDSS pipeline photometry. While more 
recent SDSS data releases (e.g. DR7; Abazajian et al. 
2009) may contain more of the optical counterparts to 
ES sources, the re-processed photometry is essential for 
accurate photometric studies and only exists for the DR2 
sample (W09). 

To be included in the ES/SDSS sample, each candi- 
date galaxy had to meet 4 criteria: 1) the ES recessional 
velocity must agree to within twice the velocity width 
(as measured by the W20 value) of the optically derived 
redshift; 2) there must be no more than 1 detectable, spa- 
tially resolved galaxy within the ES beam at the same 
redshift; 3) the candidate galaxy must not extend across 
two or more SDSS fields; and 4) all galaxies must be at 
least 1' away from any saturated foreground stars. For 
additional details on catalog matching, see W09. 

3. RESULTS 

3.1. Modeling the Colors of Hl-selected Galaxies 

Figure [T] shows the extinction corrected r — i colors of 
the galaxies in our HI selected sample as a function of 
their g — r color. The galaxies are color-coded according 
to their gas fraction. The gas fraction (f ga s) is defined 
as: 
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where Mhi is the HI mass and M* is the stellar mass de- 
termined using color-dependent mass-to-light ratios that 
have been corrected for a Kroupa initial mass function 
(Bell et al. 2003; Pizagno et al. 2005; see W09). The 
factor of 1.4 corrects for the mass in helium. We have ne- 
glected the molecular gas component, which is measured 
to be small in low-mass, LSB galaxies that dominate 
our sample (M H2 /M H i < 0.1; Schombert et al. 1990; 
Matthews et al. 2005). No molecular gas measurements 
(CO) exist for any of our sample galaxies. However, re- 
cent CO observations of nearby galaxies give some insight 
to the range molecular gas content in galaxies of similar 
morphological type. We estimate that the molecular to 
neutral gas fraction (M# 2 /M#j) may be as high at 0.5 
for some of the early-type spirals in the ES/SDSS sample 
(Leroy et al. 2009). Thus, we may be underestimating 
the gas fractions in a small number of our galaxies. 

Figure [1] shows several important trends. First, the Hl- 
selected galaxies occupy a well-defined locus in color, as 
seen in other SDSS studies (Strateva et al. 2001; Blanton 
et al. 2003a; Baldry et al. 2004; Mailer et al. 2008). Sec- 
ond, the gas fraction scales with galaxy color such that 
bluer galaxies have a larger fraction of their baryonic 
mass in gas, as has also been reported for less uniform 
samples (Kannappan 2004). Third, the scatter in color- 
color space is large for gas-rich galaxies. 
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Fig. 1. — r — i vs g — r for ES/SDSS sample galaxies. All 
photometric corrections have been applied to these data. Galaxies 
have been color coded according to their gas fractions. 



Fig. 2. — r — i vs g — r for HI selected galaxies. Galaxies 
have been color coded according to their gas fractions. Bruzual & 
Chariot (2003) population synthesis grids have been overplotted for 
various continuous SFH and metallicities. The SFH are given by 
t values for an exponential SFR starting at 12 Gyr ago. (e - */ T ). 
Metallicities are plotted in comparison to solar. 



To interpret the trends observed in Figure [TJ we com- 
pared the colors of the ES/SDSS sample galaxies to theo- 
retical predictions by running a series of Bruzual & Char- 
lot (2003) population synthesis models at various metal- 
licities and star formation histories. The models use the 
Chabrier (2003) initial mass function (IMF) to model 
the stellar populations with a range of stellar masses of 
O.6M < m < 12OM , at 6 different metallicities (0.005, 
0.02, 0.2, 0.4, 1 and 2.5 times solar) and 10 different 
star formation histories (SFH). The SFHs are modeled 
as having either (1) an exponentially decreasing star for- 
mation rate (SFR; e~'/ r ) with r values of 8, 4, 2, and 1 
Gyr (these models have mean, mass- weighted stellar ages 
of 7.4, 8.6, 10.0, and 11.0 Gyr respectively); (2) an expo- 
nentially increasing SFR with r values of 8, 4, 2, and 1 
Gyr (mean stellar ages of 4.6, 3.4, 2.0 and 1.0 Gyr respec- 
tively); (3) a continuous SFR (r = oo; mean stellar age 
of 6 Gyr) ; or (4) an instantaneous event of star formation 
(r = 0; mean stellar age of 12 Gyr). For each model we 
assume that 12 Gyr has elapsed since the first formation 
of stars (Dalcanton & Bernstein 2002; Brooks et al. 2009) 
and extract the SDSS model colors at that epoch. We 
do not include dust in any of the population models be- 
cause we have corrected for the internal extinction of the 
galaxies using the measured rotation velocities (HI line 
widths) and the prescription of Tully et al. (1998), where 
galaxies with the same rotation velocity are assumed to 
have similar internal dust extinction. The uncertainty 
in this correction is ~ 0.1 magnitudes and may increase 
the scatter of our photometry. However, many of the 
galaxies in the SDSS/ES sample have little or no dust 
correction due to their slow rotation (Tully et al. 1998; 
Dalcanton et al. 2004; see W09 for more details). 

Each Bruzual & Chariot (2003) model was run with- 
out gas recycling - the population of stars maintains the 
same metallicity throughout the star formation history. 
While this does not accurately mimic the chemical evolu- 
tion of individual galaxies (which are some amalgamation 
of ages and metallicities), it does provide a set of specific 
model comparisons by which to compare the mean prop- 
erties of observed systems. 

Figure [2] shows the resulting population grid for the 
r — i colors as a function of g — r. The "horizontal" 
lines are lines of constant metallicity and the "vertical" 



lines are lines of constant r (a proxy of mean stellar age). 
Tracks of constant age are well separated in g — r, mak- 
ing g — r a good probe of stellar age. On the other hand, 
the r — i color is only marginally sensitive to metallic- 
ity and the aforementioned degeneracy is apparent at 
all sub-solar metallicities. These grids bracket all rea- 
sonable continuous star formation histories, but bracket 
only a narrow region of g — r vs. r — i color space. Al- 
though SDSS provides photometry in 5 different photo- 
metric passbands, most of our color analysis uses only 
the g, r and i-bands. For most of our galaxies, the u and 
z bands yield photometry with large uncertainties that 
does not help us differentiate between models or con- 
strain photometric relationships. Because the it-band is 
sensitive to the star formation rate (Hopkins et al. 2003), 
we have included w-band photometry in a limited amount 
of our analysis. 

The ES/SDSS sample galaxies have been overplotted 
and color-coded according to their gas fractions, as in 
Figure 1. At the red end of the galaxy distribution, the 
colors are well explained by the model grids. The redder 
galaxies are consistent with having SFHs that have been 
declining to the present. These results are in agreement 
with previous work by Bell et al. (2003) as well as Juneau 
et al. (2005), who demonstrate that galaxies with large 
stellar masses formed most of their star in the first 3 
Gyr after their formation. The reddest galaxies also have 
colors indicative of super-solar metallicities. These red 
(g — r > 0.6), high- metallicity galaxies arc all early-type 
spiral, lenticular or elliptical galaxies. 

The data and models in Figure [5] suggest that the 
red galaxies formed almost all of their stars in the dis- 
tant past and have been undergoing little to no star for- 
mation since, in spite of containing HI at the present. 
Some mechanism is therefore responsible for the contin- 
ued presence of HI. The cannibalization of gas rich dwarfs 
or infall from "cold mode" accretion (Katz et al. 2003; 
Brooks et al. 2009) might be sufficient to explain the cur- 
rent small reservoir of gas, which is only a small fraction 
of the galaxies' baryonic mass. Regardless of the phys- 
ical process by which these galaxies acquired or main- 
tained gas, the presence of HI suggests that recent star 
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formation is likely to have taken place. Indeed, all of 
the galaxies redder than g — r > 0.5 have SDSS spectra 
that contain emission lines. In spite of the aperture bias 
in SDSS fiber spectroscopy (fibers are only 3" in diam- 
eter), these spectra can be useful for determining gross 
spectroscopic properties .Based on their location on a 
Baldwin-Phillips-Terlevich diagram (Baldwin, Phillips & 
Terlevich 1981), five of the 156 emission-line galaxies are 
consistent with hosting an AGN. All five of the galax- 
ies with AGN have PetroR90 radii (radii containing 90% 
of the Petrosian flux) larger than 30". Therefore, their 
disk emission (not the AGN) dominates the integrated 
photometry and they are included in remainder of the 
analysis. The majority of the galaxies in the ES/SDSS 
sample have emission lines that suggest that current star 
formation is underway in these galaxies. 

As we look to the bluer galaxies in Figure [H the in- 
ferred metallicities are lower and the SFHs move toward 
a constant star formation rate (r = oo). This shift is 
likely due to a decrease in star formation efficiency. We 
define star formation efficiency as the efficiency in con- 
verting gas into stars. As discussed above, we neglect the 
role of molecular gas (the state of gas from which stars 
form) , which should be irrelevant for the bluer galaxies in 
our sample. Recent studies (that included the molecular 
component) have seen a similar decrease in star forma- 
tion efficiency in nearby, Hi-dominated galaxies (Bigiel 
et al. 2008). Although we do not resolve the HI emission, 
the failure to convert large amounts of HI into stars in the 
blue galaxies constitutes a decrease in star forming effi- 
ciency. Such a decrease slows the rate of star formation, 
keeping metallicities low and gas fractions high, allowing 
stars to form at the present day. 

The bluest galaxies (g — r < 0.3) are best modeled by 
SFHs that have increasing SFRs. This area of the model 
parameter space has been ignored by previous galaxy 
studies (Kauffmann et al. 2003; Brinchmann et al. 2004; 
Salim et al. 2005; Johnson et al. 2007). Gas infall might 
serve to increase the HI surface density and subsequently 
the SFR in these galaxies. Although recent bursts have 
likely occurred in these galaxies, the already blue g — r 
colors of these do not change significantly when recent 
bursts are added to the models (see §3.3). 

Many of the bluest galaxies in the ES/SDSS sam- 
ple have colors that are not consistent with any of the 
Bruzual & Chariot models. Most of these galaxies have 
r — i colors that are too blue to be explained by any it- 
eration of the models. The LSB nature of many of these 
systems suggests that for even small SFR, emission lines 
might be able to dominate the broadband colors (Zackris- 
son et al. 2001; Magris, Binette & Bruzual 2003; Anders 
& Alvensleben 2003; Zackrisson, Bervall & Ostlin 2005; 
Zackrisson, Bergvall, & Leitet 2008). We investigate this 
possibility in detail in §3.3. 

Finally, there are a five galaxies that seem to be too red 
in r— i to be explained by any population models. Two of 
these galaxies have a saturated star nearby that is most 
likely affecting their colors. One of the galaxies has very 
large photometric uncertainties because of its low surface 
brightness and the other two have very compact cores, 
which suggests the presence of AGN. If we compute the 
colors for these two galaxies excluding the central region, 
using a shell with radial boundaries defined to be from 




Fig. 3. — r — i vs g — r for HI selected galaxies (colored sym- 
bols) plotted with a volume selected sample of galaxies from SDSS 
(contours). No internal extinction has been applied to these data 
but a reddening vector has been plotted for reference. Galaxies 
in the thesis sample have been color coded according to their gas 
fractions. Note that the HI selected sample appears to have colors 
that are shifted off of the SDSS galaxy locus. 



half of the Petrosian radius to the Petrosian radius, both 
galaxies fall back to or below the model grids when their 
central core is excluded. 

3.2. The Colors of HI vs. Optically Selected Galaxies 

Although the trends seen in Figure Q] have been noted 
in optically selected SDSS samples (e.g. Strateva et 
al. 2001; Blanton et al. 2003a), the HI selection intro- 
duces a systematic shift in the color distribution of the 
ES/SDSS sample galaxies. Figure [3] shows the r — i vs. 
g — r colors for the ES/SDSS sample plotted on top of a 
volume limited sample drawn from the SDSS DR4 sam- 
ple (Adelman-McCarthy et al. 2006). All galaxies have 
been k-corrected and corrected for Milky Way extinc- 
tion. We used the Milky Way extinction values from 
Schlegel et al. (1998) and calculated K-corrections using 
the IDL package kcorrect_v3 . 2 (Blanton et al. 2003b). 
The Blanton et al. (2003b) k-correction method uses 
empirically-derived eigenspectra to generate model spec- 
tra from SDSS photometry. All of the galaxies (both 
SDSS DR4 and ES/SDSS) were corrected using the same 
method to z=0. Most of the k-corrections were smaller 
than 0.01 magnitudes in all bands. No data in Figure [3] 
have been corrected for internal extinction. The redden- 
ing vector (derived from Schlegel et al. 1998) has been 
included to show that any correction to cither data set 
is almost perfectly aligned with the galaxy locus. 

The Hl-selected galaxies do not show the color bi- 
modality seen in previous SDSS studies (e.g. Blanton 
et al. 2003a; Baldry et al. 2004). Figure H shows the 
u — r colors as a function of stellar mass for the ES/SDSS 
galaxies (colored according to gas fraction) and the SDSS 
DR4 sample (contours). The bimodal galaxy distribution 
(blue and red sequence) is easily seen in the SDSS DR4 
galaxies, but is absent from the ES/SDSS galaxies. The 
HI selection of the ES/SDSS galaxies identifies gas-rich 
galaxies that are on the blue sequence and are currently 
in the process of star formation. Like in Figure [31 the 
ES/SDSS galaxies in Figure Q] appear to be offset from 
the SDSS DR4 sample (bluer in g — r or redder in r — i). 
Figure [4] also demonstrates the ability of an Hl-selected 
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Fig. 4. — The u — r colors as a function of stellar mass for 
the ES/SDSS galaxies (colored according to gas fraction) and the 
SDSS DR4 sample (contours). The bimodal galaxy distribution 
(blue and red sequence) is easily seen in the SDSS DR4 galaxies, 
but absent from the ES/SDSS galaxies. The HI selection of the 
ES/SDSS galaxies identifies gas-rich galaxies that are on the blue 
sequence and are currently in the process of star formation. 



sample to recover galaxies with small stellar masses; the 
ES/SDSS sample has a significant number of galaxies 
with stellar masses below 10 8 M©. 

It is clear from Figure [3] that the HI selected sources 
fall to one side of the SDSS galaxy distribution, such that 
they are either bluer in g — r or redder in r — i. This offset 
is particularly pronounced at the red end of the color 
distribution where the color dispersion is small. As noted 
in the previous section, the redder galaxies have SFHs 
that suggest that almost all star formation occurred in 
the distant past. The presence of HI in these systems 
may indicate a recent acquisition of HI and a subsequent 
ignition of star formation, as confirmed by their SDSS 
spectra. This temporary increase in SFR due to late- 
time accretion would naturally explain why the red HI 
selected galaxies are slightly bluer in g — r than the rest 
of the SDSS galaxy population (Figure [3]). 

To constrain the amount of recent star formation 
needed to produce the observed color offset between the 
ES/SDSS galaxies and the "main" SDSS sample, we re- 
ran the same series of Bruzual & Chariot (2003) popula- 
tion synthesis codes as in Figure [21 but with the addition 
of a recent burst of star formation. We generated an ini- 
tial model with g — r and r — i colors of 0.75 and 0.3 
respectively (r=2; Z=Z©). We then added a single in- 
stantaneous burst (delta-peak) of star formation to the 
continuous model and varied the time of the burst from 
2 Gyr to 100 Myr in the past. We also varied the burst 
strength from 0.1-10% of the integrated star formation. 
The resulting models were most sensitive to the time 
at which the burst was placed. The models that most 
closely matched the slight g — r bluing were bursts placed 
at 300 Myr in the past with 1% of the past integrated 
star formation occurring in the burst. 

Figure [5] shows the evolutionary track for the best 
matched model burst (300 Myr in the past) with symbols 
plotted every 50 Myr since burst. The figure demon- 
strates that there is a rapid and short-lived change in 
color that occurs in the first few Myrs and is followed 
by a slow reddening for hundreds of Myrs. Colors bluer 
than g — r < 0.4 and r — i < 0.2 are observed only for 
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Fig. 5. — The r — i vs. g — r evolution of a single Bruzual & 
Chariot model with the addition of an instantaneous burst of star 
formation that produced 1% of the total stellar mass of the galaxy 
and occurred 300 Myr in the past (solid line). The Bruzual & 
Chariot continuous star formation grid (from Figure [2) is plotted 
for reference (dashed lines). The continuous model to which the 
burst was added corresponds to a red galaxy from the ES/SDSS 
sample (t=2, Z=Zq; g — r and r — i colors of 0.75 and 0.3 respec- 
tively). Colored symbols are plotted every 50 Myr starting 50 Myr 
before the burst. Colors bluer than g — r < 0.4 and r — i < 0.2 are 
observed only for very short durations (<1 Myr). By 100 Myr the 
colors have reddened to within Ag — r=0.15 and Ar — i=0.05 of 
their initial value and continue to redden for hundreds of Myr. 



very short durations (<1 Myr). By 100 Myr the colors 
have reddened to within Ag — r=0.15 and Ar — i=0.05 of 
their initial value and continue to redden for hundreds of 
Myr. This "bottleneck" in the burst evolution suggests 
that the small color offset observed in Figure [3] could be 
sufficiently long lived to be due to a burst in the past few 
hundred Myr. 

Figure [S] shows the ES/SDSS galaxies and the model 
grids for continuous star formation with an additional 
1% burst occurring 300 Myr in the past. The red HI 
selected galaxies in the burst model have r ~ 3, com- 
pared to the r ~ 4 seen in the burst-free models. By 
comparing the difference between the grids in Figure [2] 
and [6] with the offset observed in Figure El it is clear 
that the HI selected galaxies in the burst case (Figure 
[6]) occupy the same model space (r ~ 3) as the SDSS 
"main" galaxies in the burst-free case (Figure |2j)- Unfor- 
tunately, our inability to correct for internal extinction 
in the SDSS "main" sample prohibits us from overlay- 
ing accurate grids on Figure [3l which would be highly 
instructive. However, "blinking" back and forth between 
Figures [5] and [5] confirms that a recent burst can explain 
the offset of the red galaxies. 

It is important to note that the addition of the burst 
may effect the colors of the blue galaxies as well. A single 
instantaneous burst produces extremely blue colors for 
very short times. However, in reality, the bursts are not 
instantaneous, but instead last for a finite amount of time 
with some exponential decrease (10-300 Myr depending 
on galaxy type; McQuinn et al. 2009). If a significant 
fraction of the star formation has occurred in the very 
recent past, then some of the blue colors may be due to 
recent bursts. However, rising SFRs are still required to 
explain the bulk of the galaxies with colors bluer than 
g-r <0.3. 

Although degeneracies between the time, duration and 
strength of the star formation events in the models pre- 
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Fig. 6. — r — % vs g — r for HI selected galaxies. All photometric 
corrections have been applied to these data. Galaxies have been 
color coded according to their gas fractions. Bruzual & Chariot 
(2003) population synthesis grids have been overplotted for various 
SFH and metallicities after 12 Gyrs of star formation. The SFH 
are given by r values for an exponential SFR (e~ '/ T ). These grids 
include an instantaneous burst of star formation that produced 1% 
of the stellar mass, 300 Myr in the past. 



vent an exact characterization of the bursts in the HI 
selected galaxies, we can safely claim that a new gener- 
ation of stars has been formed in the last few hundred 
Myr in the red systems (a large instantaneous burst has 
a similar effect on the integrated colors as a small expo- 
nentially decreasing burst). The bursts require no more 
than 1% of the mass to be involved as long as the burst 
was in the last 300 Myr. 

We note that the bluer g — r colors in the red ES/SDSS 
galaxies may not be solely due to main sequence stars, 
since most blue main sequence stars will have ceased their 
hydrogen burning after a few hundred Myr; the g — r 
color of a main sequence turnoff star at 300 Myr is ~0.0 
(Sarajedini et al. 2004; Covey et al. 2007). We suggest 
that another possibility for the blue g — r colors we see in 
the red HI selected galaxies could be due to an elevated 
level of BHeB stars left over from a recent burst of star 
formation. These stars are seen several hundred Myr 
after a star formation event (Dohm-Palmer et al. 2002). 

3.3. Modeling Galaxy Colors with Emission Lines 

Even with the addition of bursts, the colors of many 
ES/SDSS galaxies do not fall on the Bruzual & Chariot 
(2003) model grids. Unlike the red galaxies, which can 
be easily explained, the majority of the discrepant galax- 
ies are significantly bluer in r — i than any models can 
predict. These systems also have high gas fractions and 
most have low surface brightness. No matter what pa- 
rameters we alter in the Bruzual & Chariot models, we 
cannot produce a stellar population with very blue r — i 
colors. 

We show that the discrepant colors can be explained 
by including the emission lines that dominate the spectra 
of these outlying galaxies. The blue colors and high gas 
fractions of the discrepant galaxies suggest that they are 
actively forming stars. Because these galaxies also have 
low surface brightnesses, it is possible that emission lines 
produced by HII regions contribute significantly to their 
luminosity. HII regions produce numerous emission lines 
in the g and the r bands. Ha dominates in the r-band 
while [OIII] and H/3 dominate in the g-band. In contrast, 
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Fig. 7.— SDSS fiber spectra for 4 galaxies in the ES/SDSS 
sample with strong emission line features. All of the fibers for 
these spectra were placed on top of HII regions in the galaxies. 
The flux axis is in units of 10 — ergs/s/cm A. 



the i-band has almost no emission lines. Therefore, in an 
HII region, we would expect the r — i color to be signif- 
icantly bluer than the underlying stellar continuum and 
the g — r color to change only slightly. This suggests 
that i-band luminosity may be a better indicator of the 
underlying stellar mass than g — r. It also implies that 
the g — r color should be used to compute color depen- 
dent mass-to-luminosity ratios (the emission lines tend 
to cancel out; see below). 

We can test this idea using the SDSS spectroscopic 
fibers, which were often placed on HII regions for LSB 
galaxies. By comparing the colors of HII regions (derived 
from SDSS spectra) with the integrated colors of the 
galaxies that host them, we can quantify color changes 
produced by HII regions. 

We first selected 10 ES/SDSS galaxies with SDSS spec- 
troscopy and required that they have large Ha equivalent 
widths (> 300 A). Four of the spectra selected are shown 
in Figure [JJ All of the selected spectra have good signal- 
to-noise (> 50) and are dominated by emission line fea- 
tures. We visually inspected each galaxy to ensure that 
all of the spectra selected (and shown in Figure [JJ are 
from fibers placed on high surface brightness HII regions. 
We convolved the g, r and i SDSS filter curves with each 
spectrum and converted the computed flux density to AB 
magnitudes using the relation that a magnitude object 
has a flux density of 3631 Jy (Oke & Gunn 1983). We 
then compared the integrated colors of the 10 galaxies 
derived in W09 with the HII region colors derived from 
the spectra. 

Figure [8] shows the resulting differences between the 
HII region and the integrated color of each galaxy. The 
colors derived from fibers placed on HII regions are plot- 
ted as black symbols and the integrated galaxy colors 
derived from images (W09) are plotted as red diamonds. 
Blue lines connect a galaxy's integrated color to the color 
of the 3" diameter aperture around the HII regions. The 
HII region colors are bluer in r — i by ^0.5 m and all but 
2 become bluer in g — r by ~0.2 m , compared to the in- 
tegrated colors. This result confirms the hypothesis that 
emission lines can make the r — i color blue enough to 
explain the colors we see in the ES/SDSS sample. How- 
ever, contrary to our initial expectation, the g — r color 
does change considerably. The bluer g — r colors near HII 
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Fig. 8. — r — i vs. g — r color difference between the integrated 
colors of 10 galaxies (red diamonds) and the "fiber" magnitudes 
computed from the HII region spectra (black symbols). The blue 
lines connect each galaxy's integrated color to its HII region color. 
All of the galaxies get bluer in r — i and all but 2 get bluer in g — r 




q-r 

Fig. 9. — The r — i vs. g — r colors of the ES/SDSS galaxies with 
the addition of emission lines to an underlying stellar population 
for a single Bruzual & Chariot model galaxy (t=8 Gyr; Z=0.4 Zq). 
The colored symbols represent different star formation rates. The 
galaxy's r — i color becomes bluer as star formation is increased, 
while the g — r color goes almost unchanged. 



regions can be explained as being produced by a young 
underlying stellar population. More surprisingly, the col- 
ors near two of HII regions are redder in g — r than the 
integrated light from the galaxies. This difference may 
reflect the presence of dust localized to the HII region. 

To further test the emission line hypothesis as well as 
to confirm the reason for the g — r color change in Fig- 
ure [3 we added model emission lines to the computed 
Bruzual & Chariot model colors. Theoretical models for 
line emission in galaxies are given by Kewley et al. (2001) 
using the STARBURST99-MAPPINGS III code. The mod- 
els cover a range of metallicities, ionization states, elec- 
tron densities, and assumption about the burst lengths. 
For each model, the wavelength of every line and its en- 
ergy are provided. We convolved these energies with the 
SDSS filters and converted the computed flux densities 
into AB magnitudes. We added the emission line output 
to the Bruzual & Chariot colors using a scaling coeffi- 
cient that determines the strength of the emission line 
relative to the true stellar continuum. This coefficient is 
directly proportional to the SFR and can be used to es- 
timate the range of possible SFRs for a galaxy. The final 
models are not fully self consistent, because the ionizing 
radiation for the emission line spectrum does not come 
directly from the Bruzual & Chariot spectrum. How- 
ever, because the star formation is highly localized in 
late-type galaxies, it is possible for the ionizing radiation 
in isolated HII regions to be different than that of the 
broader component of stars 

When adding the STARBURST99-MAPPINGS III 
emission-line models to the Bruzual & Chariot outputs, 
we ensured that the adopted gas metallicity matched 
the assumed stellar metallicity. Because there are only 3 
overlapping metallicities, we limited our analysis to 0.2 
Z , 0.4 Z Q and solar metallicity populations. In reality 
most of the galaxies where emission lines dominate are 
low mass and gas rich and thus are unlikely to have high 
metallicities. Even values of 0.2 Z Q might be too high 
for these systems (Tremonti et al. 2004). 

Figures [9] and [H)] show the colors (r — i vs.g — r and 
u — r vs. g — r respectively) that result from adding 
emission lines to a single Bruzual & Chariot model (r=8 
Gyr; Z=0.4Z Q ). This particular model was chosen be- 




FlG. 10. — Theu-rvs. g — r colors of the ES/SDSS galaxies with 
the addition of emission lines to an underlying stellar population 
for a single Bruzual & Chariot model galaxy (t=8 Gyr; Z=0.4 Zq). 
The colored symbols represent different star formation rates. 



cause it lies on the red end of the blue population and 
intersects the ES/SDSS color-color locus. The colored 
symbols represent different star formation rates. In Fig- 
ure [31 the model galaxy's r — i color becomes bluer as 
star formation is increased, while the g — r color is al- 
most unchanged. A similar effect is seen in the u — r color 
in Figure [TU1 This confirms that the non-zero g — r color 
offsets near HII regions must indeed be due to stellar 
populations and or dust. 

In Figure 111) we show stellar population grids similar 
to those in Figure^ but with an 8 Myr long continuous 
burst of star formation. Black lines indicate emission 
line models with an electron density of 10 cm -3 and an 
ionization parameter of 5xl0 6 cm s _1 . Red lines have 
the same electron density but an ionization parameter 
of 8xl0 7 cm s _1 . This plot demonstrates that adding 
emission lines of varying strength can explain all of the 
discrepant blue galaxies. While 8 Myr may be short for 
some realistic starburst events, it successfully demon- 
strates the effect of current (or recent) star formation 
on the colors of galaxies. 

We found that varying the electron density and burst 
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Fig. 11. — r — i vs. g — r colors for the HI selected galaxies. 
Color coding indicates the gas fractions of the galaxies. The grids 
represent 2 different linear combinations of the STARBURST99- 
MAPPINGS III emission line models and the Bruzual & Chariot 
population synthesis models. The Bruzual & Chariot models arc 
identical to those shown in Figure [2] All emission line models use 
an 8 Myr continuous burst of star formation. Black lines indicate 
emission line models with an electron density of 10 cm -3 and an 
ionization parameter of 5xl0 6 cm s _1 . Red lines have the same 
electron density but an ionization parameter of 8xl0 7 cm s~ x . 
Solid lines have SFR densities of 1.2xl0 -3 Mq yr" 1 kpc -2 and 
3.7xl0 — 3 Mq yr — 1 kpc -2 for the two models (black and red) 
respectively. The dashed lines have SFRs 3 times larger than the 
solid lines. These models clearly describe the colors of all of the 
galaxies that are not fit by normal populations of stars. 



length (8 Myr continuous vs. instantaneous) did not 
have a major effect on the output models. All signif- 
icant changes in color are driven by the ionization pa- 
rameter and by the global normalization of the emission 
line spectrum, which is directly proportional to the SFR 
(the emission lines are generated for a specific value of 
the SFR). Solid lines in Figure fTTI have SFR densities 
of 1.2xlCT 3 Mq yr" 1 kpc" 2 and 3.7xl0" 3 M yr" 1 
kpc" 2 for the two models (black and red) respectively. 
The dashed lines have SFRs 3 times larger than the solid 
lines. For typical galaxies in our sample, we estimate 
SFRs of 0.004-0.2 M yr" 1 in the emission line dom- 
inated regions. These are reasonable SFRs for small, 
blue galaxies (van Zee 2001; Hunter & Elmcgreen 2004; 
Salzer et al. 2005). We estimated the SFR density of 
the ES/SDSS sample using the it-band luminosity as a 
proxy for SFR (Hopkins et al. 2003). We found that the 
SFR densities range from 0.0001 to 0.02 M Q yr" 1 kpc" 2 . 
This range is consistent with the star formation rate den- 
sities required to produce the blue offset in the galaxies 
in Figure [9] While the models might not be perfectly 
self-consistent as discussed above, they do demonstrate 
that strong emission lines can certainly affect the col- 
ors of galaxies. We note that the "pitched roof" shaped 
grids are a result of how the emission lines are affected 
by slight changes in the metallicity. It is clear that the 
span of these models can easily explain the galaxies that 
the stellar population models cannot and that many low 
mass, blue galaxies must have colors that are dominated 
by emission lines (e.g. Zackrisson et al. 2005). 

The emission line effect may also amplify the underly- 
ing color offset discussed in the previous section by mak- 
ing galaxies appear bluer in r — i than they actually are. 
However, the colors of the red galaxies in the offset re- 
gion are dominated by their stellar populations and not 
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Fig. 12. — Principal axis of the galaxy color-color relation af- 
ter PCA. Galaxies are color coded according to their gas fraction. 
Error bars were calculated using the uncertainties in the g, r and 
i-band photometry. 



likely affected by the emission line colors. 

3.4. The Increase in Color Dispersion 

The analysis above demonstrates that the observed col- 
ors of the ES/SDSS sample can be well explained by 
stellar population models including emission lines, and 
that the general trend in color is due to decreasing mean 
stellar age with increasing gas fraction. However, while 
this explains most of the general trends seen in the color 
distribution in Figure [TJ it does not fully explain the in- 
crease in the dispersion of color from red to blue galax- 
ies. The narrow dispersion at red colors grows to a much 
larger value before the photometric uncertainties become 
large enough to be solely responsible for the discrepancy. 
Although some of the increased dispersion is likely due to 
the onset of emission line dominated colors, this cannot 
account for the entirety of the color dispersion. 

To explore how the dispersion in color relates to other 
physical quantities, we performed a principal component 
analysis (PCA) on the r — i vs. g — r colors. The result 
aligns the principal axis (PI) with the galaxy color-color 
locus and sets the secondary axis (P2) perpendicular to 
PI. PI and P2 can be expressed in term of the g — r and 
r — i colors of the galaxies by 



PI = 0.6(r-i)+0.8(g-r) 
and 



(0.6(r 0.8(5 -r)) (2) 



P2 = O.8(r-i)-O.6(0-r)-(O.8(r-i)-O.6(s-r)), (3) 

where the angle brackets indicate the mean of the en- 
closed quantity for the sample (the resulting eigenvalues 
are 6.2 and 0.8 for the eigenvectors pointing along the 
PI and P2 axis respectively). 

Figure [12] shows the result of the PCA coordinate 
transformation. Error bars were calculated using the un- 
certainties in the <?, r and i-band photometry. The P2 
axis serves as an indicator of dispersion from the galaxy 
locus and we compare the P2 values to other physical 
parameters below. 

We quantified the dispersion in the P2 axis by calcu- 
lating the rms P2 scatter in bins of PI. Figure [TBI shows 
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Fig. 13. — The rms scatter in the P2 axis as a function of PI. 
The error bars indicate the combined photometric error of the data 
points in each bin. 



0.3 
0.2 
0.1 
0.0 




M 


-0.1 


* *** 




-0.2 
-0.3 
-0.4 




(0) 


10' 10° 10 
Stellor 


10'° 10" 1C 
MOSS (M s ) 












■ • 
















* ■ M 





24 22 20 IS 
-bond surfoce brightness (mog/orcsecond 2 ) 



1 10 100 

Stellor Surfoce Density (Ug/pc 1 ) 



Fig. 14. — P2 PCA axis as a function of stellar mass (a), velocity 
width (b), i-band surface brightness (c) and stellar surface density 
(d). Points have been color coded according to their gas fractions. 



the resulting change in color dispersion along the princi- 
pal axis. The error bars indicate the photometric uncer- 
tainty in the P2 value. As seen in Figures [T2l and [Tl blue 
galaxies have a significantly larger color dispersion than 
red galaxies. 

Changes in observable properties of the galaxies, such 
as stellar mass, rotation velocity, surface density and sur- 
face brightness might be related to a change in their color 
dispersion. As seen in previous work on the ES/SDSS 
sample, many of these physical properties are related 
to each other (Disney et al. 2008; Garcia-Appadoo et 
al. 2009) and a change in one property is likely to show 
up as a variation in the others. 

To examine how the change in color dispersion relates 
to the physical properties of the ES/SDSS galaxies, we 
plotted the P2 PCA values as a function of stellar mass, 
rotation velocity (measured from half of the 20% HI ve- 
locity width), surface brightness and stellar surface den- 
sity (calculated by measuring the average stellar mass 
density within the 90% Petrosian radius). Figure [14] 
shows that changes in the physical properties trace the 
change in dispersion very well. In addition, Figure [TJ] 
has been color coded according to gas fraction. As seen 
in previous studies, the gas fractions also correlate with 
the physical properties of the ES/SDSS galaxies (Garcia- 
Appadoo et al. 2009). We quantified the change in color 
dispersion by plotting the rms scatter in the P2 axis as 
a function of the four physical properties (Figure I15j) . 
All four properties correlate with the rms scatter in the 
P2 variable. We computed the Spearman's rank corre- 
lation coefficient for each relation to test which prop- 
erty best correlates with the color dispersion. We calcu- 
lated Spearman's p values of 0.92, 0.82, 0.89 and 0.96 for 
the stellar mass, velocity, surface brightness and surface 
density respectively. The stellar density shows the best 
correlation with color dispersion with the lowest density 
galaxies having the highest color dispersion. 

The increase in color dispersion as a function of surface 
brightness and stellar surface density may probe a change 
in the star formation process. At the low surface bright- 
nesses, the stochastic nature of star formation (Stinson 
et al. 2007), and the dominance of emission lines in low- 
mass galaxies, controls the observed colors and creates a 
larger color dispersion, while the colors of higher surface 




24 23 22 21 20 1 
-bond Surfoce Brightness (mog/orcsec 2 ) 



Stellor Surfoce Density (Ug/pc 1 ) 



Fig. 15. — The rms scatter in the P2 axis as a function of stellar 
mass (a), velocity width (b), i-band surface brightness (c) and 
stellar surface density (d). The error bars represent the combined 
photometric error in each bin. 



brightness systems are dominated by stellar populations 
(and to a smaller extent recent bursts of star formation, 
as discussed above). Because the surface brightness and 
stellar surface density are tied closely to other proper- 
ties of the ES/SDSS galaxies (Disney et al. 2008), this 
change in color dispersion may be evidence for a signifi- 
cant change in the way stars form in galaxies. 

While the stellar mass, surface brightness and stel- 
lar surface density change smoothly, the velocity width 
shows a statistically significant (4 a) jump at a velocity 
of 80 km s _1 . Previous studies have found that the ro- 
tation velocity in galaxies correlates with transitions in 
their disk properties. Dalcanton, Yoachim & Bernstein 
(2004) showed that the presence of dust lanes in a sam- 
ple of edge-on disk galaxies is related to rotation speed. 
At V c < 120 km s _1 , no dust lanes are present in their 
galaxies. Disk instabilities occurring at V c > 120 km s" 1 
act to collapse the dust to a small scale height and pro- 
duce a larger optical depth. Star formation efficiency is 
also increased above this velocity threshold, affecting the 
SFH and colors of the faster rotators (Verde et al. 2002; 
Dalcanton et al. 2004). Lee et al. (2007) find an addi- 
tional transition around 50 km s^ 1 , where the dispersion 
in Ha EW dramatically increases. This transition is in- 
terpreted as a possible change from large-scale star for- 
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mation in spiral galaxies to stochastic star formation in 
irregular systems (e.g. Stinson et al. 2007). We observe 
a transition in the color dispersion at 80 km s _1 (Figure 
[T5b) that lies in between the Dalcanton and Lee tran- 
sitions. The transition in the ES/SDSS sample may be 
a result of a similar physical process, namely the change 
from well-organized star formation in massive disks to an 
irregular stochastic process in lower mass systems 

To illustrate the differences between systems across the 
velocity transition, we randomly selected 5 gri compos- 
ite images of galaxies from the low velocity (V ro t < 80 
km s _1 ) population (top row of Figure fl"6|) and 5 images 
from the high velocity population (V ro t > 80 km s -1 ; 
bottom row of Figure ITB]) . Figure [TC] reveals a clear mor- 
phological distinction between the two populations, with 
the higher velocity systems having higher surface bright- 
nesses, earlier morphological types and a more organized 
distribution of stars (an indication of large scale star for- 
mation) . 

4. DISCUSSION 

We used the population synthesis models of Bruzual 
& Chariot (2003) to model the SFHs and metallicities of 
galaxies in the ES/SDSS sample. We found that red 
galaxies have super-solar metallicities and have SFHs 
that are consistent with them forming the bulk of their 
stars in the distant past. These red galaxies have presum- 
ably exhausted their early gas supply but have recently 
acquired gas through mergers and infall. Although the 
infalling gas is likely low metallicity and could in prin- 
ciple dilute the metallicity, the small amount of gas in 
red galaxies (low gas fractions) is not enough to make 
an observable difference in the average metallicity of the 
ES/SDSS systems. 

Bluer galaxies have lower metallicities and their mean 
stellar ages are younger (increasing r values). The gas 
fractions suggests that these blue systems have at least 
as much gas as they do stars. Previous studies have in- 
dicated that these types of galaxies are less massive and 
have disks that are stable against collapse, making their 
star formation much less efficient. They have also likely 
retained much of their initial HI and are slowly convert- 
ing it into stars as is suggested by the large (or negative) 
t values. As galaxies pass the threshold for disk stability, 
their star formation becomes sporadic and there is a clear 
dispersion in the colors due to variations in burst age, 
burst strength and superimposed emission lines. The ide- 
alistic Bruzual and Chariot models are likely only loose 
guides for these systems, as their star formation histories 
no longer follow smooth exponential functions. 

Figure fj] indicates that there are a few galaxies in the 
regime of exponentially increasing SFR (negative r val- 
ues). This is not surprising as there are few mechanisms 
that will increase the SFR of a galaxy. The likely culprit 
for these systems is the infall of gas. As gas is accreted, 
the gas densities in the galaxies will increase and local 
gravitational collapse will become more efficient. This 
will increase the SFR as well as the eventual metallicity 
of the galaxy, pushing its g — r colors to the blue (more 
recent star formation) and its r — i colors slightly to the 
red (more metals). 

We also showed that HI selected galaxies are offset from 
the SDSS galaxy locus, especially at the red end, and that 
this is likely due to bursts of star formation in the past 



few hundred Myr. The gas that induces these recent 
bursts is not primordial and is best explained by the 
accretion of gas rich dwarfs. 

The bluest galaxies in the ES/SDSS sample are not 
explained by population synthesis models alone. Their 
colors can be modeled only with the inclusion of emission 
lines. We showed that emission line spectra with reason- 
able SFRs can explain the colors of the bluest galaxies 
in our sample. 

We also showed that the distribution of galaxies at the 
red end of the color-color locus has a very small disper- 
sion that continues to rise into the blue regime. This 
change in dispersion appears to correlate with stellar 
mass, rotation velocity, surface brightness and especially 
stellar surface density. It is possible that the change in 
dispersion is evidence for a significant change in the way 
stars form in galaxies at a given mass scale. Massive 
galaxies have unstable disks and efficiently convert most 
of their primordial reservoirs of gas into stars in the first 
few Gyrs after their formation. These galaxies are bulge 
dominated and have redder, older and more massive pop- 
ulations of stars. We do see a sharp transition in the color 
distribution at a rotation velocity of 80 km s^ 1 . This 
transition may be related to previous results that have 
found sharp transitions in galaxy properties as a func- 
tion of rotational velocity (Dalcanton et al. 2004; Lee et 
al. 2007). However, further investigations are required to 
explore the various velocity transitions as a function of 
galaxy properties. 

It may also be the case that the increase in disper- 
sion is nothing more than a selection effect related to the 
surface brightness. As the surface brightness of a sys- 
tem decreases, the effect of a single burst of stars on the 
color becomes increasingly large, possibly explaining the 
increase in scatter. This effect might explain the onset 
of emission line dominated colors in the bluest galaxies 
as they can influence of the colors of galaxies devoid of 
massive stellar populations. 

We note that the change in color dispersion is not eas- 
ily seen in the volume selected SDSS data. Applying the 
HI selection identifies a low dispersion subsample of the 
SDSS galaxies, most notably at the red end. As men- 
tioned above, the red HI selected galaxies appear to be 
bluer in g — r than the "main" SDSS sample. These two 
features are likely related. It is possible that the large 
distribution in color at the red end of the SDSS main 
sample is due to the diversity of gas content. Because 
most of the red galaxies in SDSS exhausted their orig- 
inal supply of gas long ago, this dispersion may be an 
indication of the spread in time since the last major gas 
infall. We leave further discussion of the "color offset" 
to future study. 
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Fig. 16. — 5 randomly selected gri composite images of galaxies from the slowly rotating regime (Wot < 80 km s 1 ; top row) and 5 
composite images from the faster regime (V ro t > 80 km 9 ; bottom row). [A higher resolution version of this figure is available from ADS 
and the Astronomical Journal.] 



discussions while completing this project. We also thank 
John Bochanski for technical assistance and Adam Bur- 
gasser for his financial support and insistence that this 
project reach completion. 

Funding for the Sloan Digital Sky Survey (SDSS) and 
SDSS-II has been provided by the Alfred P. Sloan Foun- 
dation, the Participating Institutions, the National Sci- 
ence Foundation, the U.S. Department of Energy, the 
National Aeronautics and Space Administration, the 
Japanese Monbukagakusho, and the Max Planck Society, 
and the Higher Education Funding Council for England. 
The SDSS Web site is |http: / /www.sdss.org/| 

The SDSS is managed by the Astrophysical Research 
Consortium (ARC) for the Participating Institutions. 
The Participating Institutions are the American Mu- 



Abazajian, K., et al. 2004, AJ, 128, 502 
— . 2009, ApJS, in press 

Adelman-McCarthy, J. K., et al. 2006, ApJS, 162, 38 

Anders, P., & Fritze-v. Alvensleben, U. 2003, A&A, 401, 1063 

Balcells, M., & Peletier, R. F. 1994, AJ, 107, 135 

Baldry, I. K., Glazebrook, K., Brinkmann, J., Ivezic, Z., Lupton, 

R. H., Nichol, R. C, & Szalay, A. S. 2004, ApJ, 600, 681 
Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 
Ball, N. M., Loveday, J., Brunner, R. J., Baldry, I. K., & 

Brinkmann, J. 2006, MNRAS, 373, 845 
Barnes, D. G., et al. 2001, MNRAS, 322, 486 
Bell, E. F., & de Jong, R. S. 2000, MNRAS, 312, 497 
— . 2001, ApJ, 550, 212 

Bell, E. F., Mcintosh, D. H., Katz, N., & Weinberg, M. D. 2003, 
ApJS, 149, 289 

Bigiel, F., Leroy, A., Walter, F., Brinks, E., de Blok, W. J. G., 
Madore, B., & Thornley, M. D. 2008, AJ, 136, 2846 

Blanton, M. R., Brinkmann, J., Csabai, I., Doi, M., Eisenstein, D., 
Fukugita, M., Gunn, J. E., Hogg, D. W., & Schlegel, D. J. 2003a, 
AJ, 125, 2348 

Blanton, M. R., Eisenstein, D., Hogg, D. W., Schlegel, D. J., & 

Brinkmann, J. 2005, ApJ, 629, 143 
Blanton, M. R., et al. 2001, AJ, 121, 2358 
— . 2003b, ApJ, 594, 186 

Brinchmann, J., Chariot, S., White, S. D. M., Tremonti, C, 
Kauffmann, G., Heckman, T., & Brinkmann, J. 2004, MNRAS, 
351, 1151 

Brooks, A. M., Governato, F., Quinn, T., Brook, C. B., & Wadsley, 

J. 2009, ApJ, 694, 396 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Chabrier, G. 2003, PASP, 115, 763 



seum of Natural History, Astrophysical Institute Pots- 
dam, University of Basel, University of Cambridge, Case 
Western Reserve University, The University of Chicago, 
Drexel University, Fermilab, the Institute for Advanced 
Study, the Japan Participation Group, The Johns Hop- 
kins University, the Joint Institute for Nuclear Astro- 
physics, the Kavli Institute for Particle Astrophysics 
and Cosmology, the Korean Scientist Group, the Chi- 
nese Academy of Sciences (LAMOST), Los Alamos Na- 
tional Laboratory, the Max-Planck-Institute for Astron- 
omy (MPIA), the Max-Planck-Institute for Astrophysics 
(MPA), New Mexico State University, Ohio State Univer- 
sity, University of Pittsburgh, University of Portsmouth, 
Princeton University, the United States Naval Observa- 
tory, and the University of Washington. 



Covey, K. R., et al. 2007, AJ, 134, 2398 
Dalcanton, J. J., & Bernstein, R. A. 2002, AJ, 124, 1328 
Dalcanton, J. J., Yoachim, P., & Bernstein, R. A. 2004, ApJ, 608, 
189 

de Jong, R. S. 1996, Journal of Astronomical Data, 2, 1 

Disney, M. J., Romano, J. D., Garcia-Appadoo, D. A., West, A. A., 

Dalcanton, J. J., & Cortese, L. 2008, Nature, 455, 1082 
Dohm-Palmer, R. C, Skillman, E. D., Mateo, M., Saha, A., 

Dolphin, A., Tolstoy, E., Gallagher, J. S., & Cole, A. A. 2002, 

AJ, 123, 813 
Doyle, M. T., et al. 2005, MNRAS, 361, 34 

Driver, S. P., Allen, P. D., Graham, A. W., Cameron, E., Liske, 

J., Ellis, S. C, Cross, N. J. G., De Propris, R., Phillipps, S., & 

Couch, W. J. 2006, MNRAS, 368, 414 
Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shimasaku, K., 

& Schneider, D. P. 1996, AJ, 111, 1748 
Galaz, G., Dalcanton, J. J., Infante, L., & Treister, E. 2002, AJ, 

124, 1360 

Garcia-Appadoo, D. A., West, A. A., Dalcanton, J. J., Cortese, L., 

k. Disney, M. J. 2009, MNRAS, 394, 340 
Gavazzi, G., & Boselli, A. 1996, Astrophysical Letters 

Communications, 35, 1 
Geha, M., Blanton, M. R., Masjedi, M., & West, A. A. 2006, ApJ, 

653, 240 

Governato, F., Willman, B., Mayer, L., Brooks, A., Stinson, G., 
Valenzuela, O., Wadsley, J., & Quinn, T. 2007, MNRAS, 374, 
1479 

Gunn, J. E., et al. 1998, AJ, 116, 3040 
— . 2006, AJ, 131, 2332 

Hambly, N. C, Davenhall, A. C, Irwin, M. J., & MacGillivray, 
H. T. 2001a, MNRAS, 326, 1315 



12 



West et al. 



Hambly, N. C, Irwin, M. J., & MacGillivray, H. T. 2001b, 
MNRAS, 326, 1295 

Hambly, N. C, MacGillivray, H. T., Read, M. A., Tritton, S. B., 
Thomson, E. B., Kelly, B. D., Morgan, D. H., Smith, R. E., 
Driver, S. P., Williamson, J., Parker, Q. A., Hawkins, M. R. S., 
Williams, P. M., & Lawrence, A. 2001c, MNRAS, 326, 1279 

Hogg, D. W., Finkbeiner, D. P., Schlegel, D. J., & Gunn, J. E. 
2001, AJ, 122, 2129 

Hopkins, A. M., Miller, C. J., Nichol, R. C, Connolly, A. J., 
Bcrnardi, M., Gomez, P. L., Goto, T., Tremonti, C. A., 
Brinkmann, J., Ivezic, Z., & Lamb, D. Q. 2003, ApJ, 599, 971 

Hunter, D. A., Elmegreen, B. G., & Baker, A. L. 1998, ApJ, 493, 
595 

Ivezic, Z., et al. 2004, Astronomischc Nachrichtcn, 325, 583 
Jarrett, T. H., Chester, T., Cutri, R., Schneider, S., Skrutskie, M., 

& Huchra, J. P. 2000, AJ, 119, 2498 
Johnson, B. D., et al. 2007, ApJS, 173, 392 
Juneau, S., ct al. 2005, ApJ, 619, L135 
Kannappan, S. J. 2004, ApJ, 611, L89 

Katz, N., Keres, D., Dave, R., & Weinberg, D. H. 2003, in ASSL 
Vol. 281: The IGM/Galaxy Connection. The Distribution of 
Baryons at z=0, 185 — h 

Kauffmann, G., & others. 2003, MNRAS, 341, 33 

Kauffmann, G., White, S. D. M., Heckman, T. M., Menard, B., 
Brinchmann, J., Chariot, S., Tremonti, C, & Brinkmann, J. 
2004, MNRAS, 353, 713 

Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & 
Trevena, J. 2001, ApJ, 556, 121 

Lee, J. C, Kcnnicutt, R. C, Funes, Jose G., S. J., Sakai, S., & 
Akiyama, S. 2007, ApJ, 671, L113 

Leroy, A. K., Walter, F., Bigiel, F., Usero, A., Weiss, A., Brinks, 
E., de Blok, W. J. G., Kennicutt, R. C, Schuster, K.-F., Kramer, 
C, Wiesemeyer, H. W., & Roussel, H. 2009, AJ, 137, 4670 

Lupton, R. H., Ivezic, Z., Gunn, J. E., Knapp, G., Strauss, 
M. A., & Yasuda, N. 2002, in Presented at the Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference, 
Vol. 4836, Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series, ed. J. A. Tyson & S. Wolff, 350-356 

MacArthur, L. A., Courteau, S., & Holtzman, J. A. 2003, ApJ, 582, 
689 

Magris C, G., Binette, L., & Bruzual A., G. 2003, ApJS, 149, 313 
Mailer, A. H., Berlind, A. A., Blanton, M. R., & Hogg, D. W. 2008, 
ArXiv c-prints 

Matthews, L. D., Gao, Y., Uson, J. M., & Combes, F. 2005, AJ, 
129, 1849 

McQuinn, K. B. W., Skillman, E. D., Cannon, J. M., Dalcanton, 

J. J., Dolphin, A., Stark, D., & Weisz, D. 2009, ApJ, 695, 561 
Meyer, M. J., et al. 2004, MNRAS, 350, 1195 
Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 

Pier, J. R., Munn, J. A., Hindsley, R. B., Hennessy, G. S., Kent, 
S. M., Lupton, R. H., & Ivezic, Z. 2003, AJ, 125, 1559 



Pizagno, J., Prada, F., Weinberg, D. H., Rix, H.-W., Harbeck, D., 
Grcbel, E. K., Bell, E. F., Brinkmann, J., Holtzman, J., & West, 
A. 2005, ApJ, 633, 844 

Roberts, M. S., & Haynes, M. P. 1994, ARA&A, 32, 115 

Rosenberg, J. L., Schneider, S. E., & Posson-Brown, J. 2005, AJ, 
129, 1311 

Salim, S., et al. 2005, ApJ, 619, L39 

Salzer, J. J., Jangren, A., Gronwall, C, Werk, J. K., Chomiuk, 
L. B., Caperton, K. A., Melbourne, J., & McKinstry, K. 2005, 
AJ, 130, 2584 

Sarajedini, A., Brandt, K., Grocholski, A. J., & Tiede, G. P. 2004, 
AJ, 127, 991 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Schombert, J. M., Bothun, G. D., Impey, C. D., & Mundy, L. G. 

1990, AJ, 100, 1523 
Searle, L., Sargent, W. L. W., & Bagnuolo, W. G. 1973, ApJ, 179, 

427 

Skibba, R. A., Bamford, S. P., Nichol, R. C, Lintott, C. J., 
Andrcescu, D., Edmondson, E. M., Murray, P., Raddick, M. J., 
Schawinski, K., Slosar, A., Szalay, A. S., Thomas, D., & 
Vandcnbcrg, J. 2008, ArXiv e-prints 

Skrutskie, M. F., et al. 2006, AJ, 131, 1163 

Smith, J. A., et al. 2002, AJ, 123, 2121 

Stinson, G. S., Dalcanton, J. J., Quinn, T., Kaufmann, T., & 

Wadsley, J. 2007, ApJ, 667, 170 
Stoughton, C, et al. 2002, AJ, 123, 485 
Strateva, I., et al. 2001, AJ, 122, 1861 
Tinsley, B. M. 1972, A&A, 20, 383 
Tinsley, B. M., & Gunn, J. E. 1976, ApJ, 203, 52 
Tremonti, C. A., Heckman, T. M., Kauffmann, G., Brinchmann, J., 

Chariot, S., White, S. D. M., Seibert, M., Peng, E. W., Schlegel, 

D. J., Uomoto, A., Fukugita, M., & Brinkmann, J. 2004, ApJ, 

613, 898 

Tully, R. B., Pierce, M. J., Huang, J.-S., Saunders, W., Verheijen, 

M. A. W., & Witchalls, P. L. 1998, AJ, 115, 2264 
van Zee, L. 2001, AJ, 121, 2003 

Verde, L., Oh, S. P., & Jimenez, R. 2002, MNRAS, 336, 541 
West, A. A., Garcia- Appadoo, D. A., Dalcanton, J. J., Disney, 

M. J., Rockosi, C. M., Ivezic, Z., Bentz, M. C, & Brinkmann, J. 

2009, AJ, submitted 
Wong, O. I., ct al. 2006, MNRAS, 371, 1855 
York, D. G., ct al. 2000, AJ, 120, 1579 

Zackrisson, E., Bergvall, N., & Ostlin, G. 2005, A&A, 435, 29 
Zackrisson, E., Bergvall, N., & Leitet, E. 2008, ApJ, 676, L9 
Zackrisson, E., Bergvall, N., Olofsson, K., & Siebert, A. 2001, A&A, 
375, 814 

Zwaan, M. A., et al. 2004, MNRAS, 350, 1210 



